Draft version February 1, 2008 

Preprint typeset using LATpjX style emulatcapj v. 14/09/00 



WEAK LENSING STUDY OF GALAXY BIASING 

Henk Hoekstra 1,2,3 , Ludovic van Waerbeke 4,1 , Michael D. Gladders 2 ' 3 ' 5 , Yannick 



(N 
O 
O 
(N 

C 



> 

m 
O 

O 
(N 
O 

^ 

Of 

6 

CO 



X 



Mellier 4 ' 6 ' 3 , and H.K.C. Yee 2 ' 3 

Draft version February 1, 2008 

ABSTRACT 

We combine weak lensing measurements from the Red-Sequence Cluster Survey (RCS) and the 
VIRMOS-DESCART survey, and present the first direct measurements of the bias parameter b and 
the galaxy-mass cross-correlation coefficient r on scales ranging from 0.2 to 9.3 Mpc (which corre- 
spond to aperture radii of 1.5' — 45') at a lens redshift z ~ 0.35. We find strong evidence that both b and 
r change with scale for our sample of lens galaxies (19.5 < Rc < 21), which have luminosities around L„. 



For the currently favored cosmology (fi m = 0.3, O a = 0.7), we find b — 0.71 



+0.06 

-0.04 



confidence) on a 



scale of 1 — 2/i^q 1 Mpc, increasing to ~ 1 on larger scales. The value of r has only minimal dependence 
on the assumed cosmology. The variation of r with scale is very similar to that of 6, and reaches a 
minimum value of r ~ 0.571qq7 (at lh^ Mpc; 68% confidence). This suggests significant stochastic 
biasing and/or non-linear biasing. On scales larger than ~ 4/iljj, 1 Mpc the value of r is consistent with a 
value of r = 1. In addition we use RCS data alone to measure the ratio b/r on scale ranging from 0.15 
to 12.5 Mpc (1' — 60') and find that the ratio varies somewhat with scale. We obtain an average 
value of b/r — 1.090 ± 0.035, in good agreement with previous estimates. A (future) careful comparison 
of our results with models of galaxy formation can provide unique constraints, as r is linked intimately 
to the details of galaxy formation. 

Subject headings: cosmology: observations — dark matter — gravitational lensing 



1. INTRODUCTION 

The growth of structures in the universe via gravita- 
tional instability is an important ingredient in our under- 
standing of galaxy formation. However, the connection to 
observations is not straightforward, as we need to under- 
stand the relation between the dark matter distribution 
and the galaxies themselves. Galaxy formation is a com- 
plex process, and it is not guaranteed a priori that this 
relation, referred to as galaxy biasing, is a simple one. 
The bias might be non-linear, scale dependent or stochas- 
tic. In the simplest case, linear deterministic biasing, the 
relation between the dark matter and the galaxies can be 
characterized by a single number b (e.g., Kaiser 1987). 

Most observational constraints of biasing come from dy- 
namical studies (see Strauss & Willick 1995) which probe 
relatively large scales (10/igQ 1 Mpc or more). Recent esti- 
mates on these scales suggest values of b ~ 1 for galaxies 
(e.g., Peacock et al 2001; Verde et al. 2001). On smaller 
scales some constraints come from measurements of the 
galaxy two-point correlation function, which is compared 
to the (dark) matter correlation function computed from 
numerical simulations. 
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These studies indicate that the bias parameter b ~ 0.7 
on scales less than ~ 2ft,^~ Mpc (f2 m = 0.3, J7a = 0.7). On 
larger scales b increases to a value close to unity (Jenkins 
et al. 1998). 

Although this procedure provides useful information 
about the bias parameter 6 as a function of scale, it does 
rely on the assumptions made in the numerical simula- 
tions. In addition, it cannot be used to examine how tight 
the correlation between the matter and light distribution 
is. To do so, we need to measure the galaxy-mass cross- 
correlation coefficient r, which is a measure of the amount 
of stochastic and non-linear biasing (e.g., Pen 1998; Dekel 
& Lahav 1999; Somerville et al. 2001). 

The coherent distortions in the shapes of distant galax- 
ies caused by weak gravitational lensing provide a unique 
way to study the dark matter distribution in the universe. 
Weak lensing probes the dark matter distribution directly, 
regardless of the light distribution. In addition, it pro- 
vides measurements on scales from the quasi-linear to the 
non-linear regime, where comparisons between observa- 
tions and predictions are still limited. Much progress has 
been made in recent years, and the latest results give ac- 
curate joint constraints on the matter density fi m , and the 
normalisation of the power spectrum as (e.g., Bacon et al. 
2002; Hoekstra et al. 2002b; Refregier et al. 2002; van 
Waerbeke et al. 2002). 

Although redshift surveys can be used to determine 
the relative values of b and r for different galaxy types 
(Tegmark & Bromley 1999; Blanton 2000), weak gravita- 
tional lensing provides the only direct way to measure the 
galaxy-mass cross-correlation function (e.g., Fischer et al. 
2000; Wilson, Kaiser, & Luppino 2001; McKay et al. 2001; 
Hoekstra, Yee & Gladders 2001b). Fischer et al. (2000) 
used the SDSS commissioning data to measure the galaxy- 
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mass correlation function and their results suggested an 
average value of b/r ~ 1 on submegaparsec scales. This 
approach has been explored by Guzik & Seljak (2001) who 
used semi-analytic models of galaxy formation combined 
with N-body simulations. Their results suggest that the 
cross-correlation coefficient is close to unity. The galaxies 
used in their analysis, however, are massive (and conse- 
quently luminous) because of the limited mass resolution 
of the simulations. 

In this Paper we use the method proposed by Schnei- 
der (1998) and Van Waerbeke (1998), which allows us to 
measure the biasing parameters directly as a function of 
scale. Hoekstra, Yee, & Gladders (2001a) applied this 
technique to 16 deg 2 of i?c~band imaging data from the 
Red-Sequence Cluster Survey (RCS; e.g., Yee & Gladders 
2001), and measured the ratio b/r as a function of scale. 
They found an average value of b/r = 1.05±g;^ for a 
ACDM cosmology on scales ranging from IbOh^ kpc to 
3/150 Mpc. 

However, an accurate measurement of the mass auto- 
correlation function is required to constrain b and r sepa- 
rately. Furthermore the analysis is facilitated if the mea- 
surements of both the galaxy auto-correlation function and 
the mass auto-correlation function probe the power spec- 
trum at the same redshift. To this end, the measurements 
presented by van Waerbeke et al. (2002), based on deep 
imaging data from the VIRMOS-DESCART survey, are 
ideal. In this Paper we combine the results of the RCS 
and the VIRMOS-DESCART survey, allowing us for the 
first time to measure both b and r as a function of scale. 

The structure of the paper is as follows. In §2 we dis- 
cuss the data sets we use for our analysis. In §3 we define 
the bias parameters, and show how they can be related to 
observed correlation functions. The relevant correlation 
functions are discussed in §4. The actual measurements 
are discussed in §5. The inferred bias parameters are pre- 
sented in §6. In the appendix we demonstrate that the 
method used in this paper gives accurate results if the 
bias paramters vary with scale and redshift. 



2. DATA 

The Red-Sequence Cluster Survey 1 is a galaxy cluster 
survey designed to provide a large number of clusters with 
0.1 < z < 1.4 (e.g., Yee & Gladders 2001). Here we use 
the Rc band data from the Northern half of the survey, 
which consists of 10 widely separated patches on the sky, 
observed with the CFH12k camera on the CFHT. The to- 
tal area observed is 45.5 deg 2 , but because of masking the 
effective area is somewhat smaller, with a total of 42 deg 2 
used in the lensing analysis. The data and weak lensing 
analysis are described in detail in Hoekstra et al. (2002a). 

The VIRMOS-DESCART survey 2 consists of four 
patches in 4 colors, also observed with the CFHT12k cam- 
era. The final survey will cover 16 deg 2 , but here we use 
8.5 deg 2 of /-band data described in Van Waerbeke et al. 
(2001). 

These data have been used elsewhere to derive joint con- 
straints on Q, m and ag using the weak lensing signal caused 

1 http: //www. astro .utoronto . ca/~gladders/RCS 

2 http : //www . astr sp-mrs . f r/ virmos/ and http : //terapix .iap.fr/ 



by large scale structure (Hoekstra et al. 2002a; Hoek- 
stra et al. 2002b; Pen et al. 2002; van Waerbeke et al. 
2001; van Waerbeke et al. 2002). These papers describe 
in detail the object detection, the shape measurements 
and the corrections for the various observational distor- 
tions (PSF anisotropy, seeing, and camera shear). These 
studies have shown that the contamination of the lensing 
signal by residual systematics is small. 

The accuracy of the measurement of the galaxy auto- 
correlation function increases with survey area, and hence 
the RCS data are best suited for this measurement. Cur- 
rently, the VIRMOS-DESCART survey provides the most 
accurate measurement of the matter auto-correlation func- 
tion (in particular on scales less than 10 arcminutes, 
where systematics are negligible), although the RCS data 
also provide useful constraints on cosmological parameters 
(Hoekstra et al. 2002b). On scales less than 10 arcminutes, 
the accuracy of the RCS measurements is likely limited by 
intrinsic alignments of the sources, whereas the VIRMOS- 
DESCART measurements do not suffer from this. For the 
galaxy-mass cross-correlation function we use the results 
from the RCS, because the lenses used for the determina- 
tion of the galaxy auto-correlation function were selected 
from this survey. 

Currently, we do not have redshift information for the 
RCS galaxies, and we therefore select a sample of lenses on 
the basis of their apparent i?c~band magnitude: we define 
a sample of lenses from the RCS with 19.5 < Rc < 21, 
which yields a total of <~ 1.2 x 10 5 lenses. For the mea- 
surement of the galaxy-mass cross-correlation function we 
use source galaxies with 21.5 < Rc < 24 which yields 
<~ 1.5 x 10 6 sources. The weak lensing analysis of the 
VIRMOS-DESCART data uses - 5 x 10 5 galaxies with 
Iab < 24.5. 

In order to interpret the measurements, we have to know 
the redshift distributions of the lens galaxies and the two 
populations of source galaxies. 

The CNOC2 Field Galaxy Redshift Survey (e.g., Yee et 
al. 2000) has determined the redshift distribution of galax- 
ies down to a nominal limit of Rc — 21.5. Hence, the re- 
sults from CNOC2 provide an excellent measure of the red- 
shift distribution of our sample of lenses (19.5 < Rc < 21), 
for which we obtain a median redshift of z = 0.35. 

The source galaxies from the RCS and the VIRMOS- 
DESCART surveys have different redshift distributions, 
and these galaxies are generally too faint for spectro- 
scopic surveys. Fortunately, photometric redshifts work 
well, as demonstrated by Hockstra et al. (2000a). We 
use photometric redshift distributions derived from both 
Hubble Deep Fields (Fernandez-Soto et al. 1999), and 
multi-colour observations done with the Very Large Tele- 
scope (Van Waerbeke et al. 2002). We find a median 
redshift of z = 0.53 for the RCS source galaxies (with 
21.5 < R c < 24). The galaxies from the VIRMOS- 
DESCART survey are fainter, and the selection of galaxies 
with Iab < 24.5 gives a median redshift of z ~ 0.9. Al- 
though the source redshift distributions are quite different, 
we will show below that the two surveys are actually well 
matched, as both probe the power spectrum at z ~ 0.35 
(see Fig. 1). 
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3. METHOD 

To study the galaxy biasing we use a combination of the 
galaxy and mass auto-correlation functions, as well as the 
cross-correlation function. The statistic we use to present 
the results is the "aperture mass" M ap , which is described 
in detail in Schneider et al. (1998). It is defined as (Kaiser 
et al. 1994) 



M ap {9) = / d^u^Mt). 



(1) 



Provided U(4>) is a compensated filter, i.e., J d<fr(f>U ((f) = 
0, with U((j>) = for 4> > 9, the aperture mass can be 
expressed in term of the observable tangential shear 7 t us- 
ing a different filter function Q{(j>) (which is a function of 

M ap {9) = f d 2 ^Q{4>) lt {4). (2) 
Jo 

We use the filter function suggested by Schneider et al. 
(1998) 



U(4>) 



ap 



with the corresponding Q{4>) 



Q{4>) 



(3) 



(4) 



A common definition of the "bias" parameter is the ratio 
of the variances of the galaxy and dark matter densities, 
which is the definition we will use here. In the case of de- 
terministic, linear biasing the galaxy density contrast 5 g is 
simply related to the mass density contrast S as 5 g = b5 
(Kaiser 1987), and the ratio of the variances is the only 
relevant parameter. The galaxy number density contrast 
An g is then given by 

An g (0) = N{9) - - =bjdw Pf (w)6(f K (w)9; w), (5) 

where TV is the average number density of lens galaxies, w 
is the comoving distance, fx{w) is the comoving angular 
diameter distance, and where pf(w)dw corresponds to the 
rcdshift distribution of lens galaxies. 

Then the "aperture count" Af is given by (Schneider 1998) 



Af(8 ap ) = J d 2 0C/(0)An g (0). 



(6) 



With our choice of the filter function U (4>) , we can write 
the auto-correlation of Af as (Schneider 1998; Van Waer- 
beke 1998; Hoekstra et al. 2001) 

(Af 2 (9 ap )) = 27Tb 2 J dwh^w; 8 ap ) (7) 
where hi(w;9 ap ) is defined as 

' Pfi w ) X " 



hi(w;0 ap ) 



Pfilter(w;0ap), (8) 



and where the "filtered" power spectrum Pfuter(w,9 ap ) is 
given by (e.g., Schneider et al. 1998) 



ffiiter(w;i 



= / dl IP 



3d 



I 



;w)J 2 (W ap ). (9) 



Jk{w)'' 

Here P^ is the time-evolving 3-D power spectrum. As has 
been shown by Jain & Seljak (1997), it is important to use 
the non-linear power spectrum in the calculations, and in 
the following we use the results from Peacock & Dodds 
(1996). The filter function J(rj) is given by 



J(v) = ^5 MO) 



(10) 



where Ji{x) is the fourth order Bessel function of the first 
kind. 

The power spectrum P^d contains a wealth of informa- 
tion about the cosmological parameters. A measurement 
of (AT 2 ) could be a powerful tool, provided the value of 
b is known. Unfortunately, the latter is not true. This 
is the reason why weak lensing by large scale structure 
("cosmic shear") has become an important cosmological 
tool: it probes the (dark) matter power spectrum directly, 
without having to rely on the light distribution. 
The matter auto-correlation function (M ap ) 2 is related to 
the power spectrum as (e.g., Schneider et al. 1998) 



«(M> = y 



n 2 



dwh 2 (w;9 lip ), (11) 



where fl m is the density parameter, and h,2(w; ap ) is given 

by 



a(w) 



Pfilter(lO;0ap), 



(12) 



where a(w) is the cosmic expansion factor. The function 
g(w) is given by 

#)=/ WH ^W)^, d3) 
Jw Jk(w') 

and depends on the redshift distribution of the (back- 
ground) sources pb(w)dw. It measures the "lensing 
strength" of a lens at a distance w. If the lens is close to the 
sources, the lensing signal decreases, whereas a large dis- 
tance between the lens and the sources results in a larger 
signal. Hence, g(w) declines with increasing w, reaching a 
value of 1 for w = 0, and a value of for lenses behind the 
sources. 

We use Eqns. 7 and 11 to relate the bias parameter b to 
the cosmology and measurements, and obtain 



b 2 = 



Ho 

c 



J dwh 2 (w;9 ap ) 
J dwhi(w;9 ap ) 



nt x 



fi(o &p ,ii m ,si\) x n 2 r< 



(M 2 p (e ap )) 

<A^ ap )> 



(M? 



ap(#ap)) 



(14) 



The value of /i depends on the assumed cosmological 
model and the redshift distributions of the lenses and the 
sources. Hence, for a given cosmology, the bias parameter 
b can be determined from the observed ratio of the galaxy 
and matter auto-correlation functions. Calculations of /i 
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as a function of aperture size show that it depends mini- 
mally on the aperture size, and the adopted power spec- 
trum. 

The bias relation is likely to be more complicated than 
the simple case of linear deterministic biasing: the actual 
relation depends on the process of galaxy formation, and 
might be stochastic, non-linear or both. Hence, there is 
no reason that b is constant with scale, but as we demon- 
strate in Appendix A, we can still use Eqn. 12 as long as b 
varies slowly with scale. Another complication is that b is 
likely to depend on redshift, and conseqently the derived 
value for b is a redshift averaged value for the lenses in our 
sample. 

If the bias parameter depends on redshift and scale, 
it is important that both the galaxy and matter auto- 
correlation functions probe the power spectrum at the 
same "effective" redshift. In addition, the interpretation 
is facilitated if the measurements probe a relatively small 
range in redshift. To examine this in more detail, it is 
useful to plot the hi(w;9 ap ) and h2(w;6 ap ) as a function 
of w. 

The results are presented in Figure 1. The integrands 
have been normalized to a peak value of unity. The thick 
solid line shows hi(w; 9 ap ) as a function of w for a fiducial 
$ap = 5'. We used the observed redshift distribution of 
galaxies with 19.5 < R c < 21 from the CNOC2 survey, 
which gives rise to the wiggles in the function. This result 
indicates that we probe the power spectrum at an effective 
redshift ~ 0.35. The thick dashed line shows h,2(w;0 ap ), 
for the VIRMOS-DESCART survey. This function is 
rather broad, which reflects the fact that weak lensing 
probes the power spectrum over a relatively large redshift 
range. However, comparison with the integrand for the 
galaxy auto-correlation function shows that both peak at 
the same redshift. Hence, we can use the ratio of (Af 2 ) 
and (M ap ) measured at the same aperture size 9 ap . 

The lensing is maximal when the the lens is half-way 
between the source and the observer. Hence, the weight 
function h,2(w) reaches its peak value approximately half- 
way between the observer and the source. Compared to 
the background galaxies from the VIRMOS-DESCART 
survey, the RCS sources are at lower redshift, and con- 
sequently the latter will probe the power spectrum at a 
lower redshift. This is indicated by the thin dashed line 
in Figure 1, which shows h2(w;6 ap ) using the RCS back- 
ground galaxy redshift distribution. It is a poor match 
to the galaxy auto-correlation function. We therefore 
use the mass auto-correlation function measured from the 
VIRMOS-DESCART survey. 

So far, we have only used the auto-correlation func- 
tions. We can, however, also measure the galaxy-mass 
cross-correlation function (AfM ap ), which can be used to 
quantify how well the mass distribution correlates with the 
light distribution. 

The galaxy-mass cross-correlation function (J\fM ap ) is 
related to the power spectrum as (Schneider 1998; Van 
Waerbeke 1998; Hoekstra et al. 2001) 



{M ap {9 ap )M{9 ap )) = 3tt i^j Q m br J dwh 3 (w;9 ap ), 

(15) 

where h 3 (w;9 ap ) is defined as 



redshift 

0.00 0.27 0.58 0.96 1.46 




Fig. 1. — The different integrands of the integrals in Eqns. 14 
and 17 for a fiducial aperture size 6> ap = 5'. The integrands have 
been normalized to have the same maximum value. The thick solid 
line shows h\(w) as a function of w. The upper axis indicates the 
corresponding redshift. We used the observed redshift distribution 
of galaxies with 19.5 < Rq < 21 from the CNOC2 survey, which 
gives rise to the wiggles in the function. This result indicates that 
we probe the power spectrum (and the bias parameters) at an ef- 
fective redshift ~ 0.35. The thick dashed line shows h2(w), for the 
VIRMOS-DESCART survey. This function is rather broad, which 
reflects the fact that weak lensing probes the power spectrum over 
a relatively large redshift range. However, comparison with the in- 
tegrand for the galaxy auto-correlation function shows that both 
peak at the same redshift. Hence, we can use the ratio of (A/ -2 ) and 
(Aff p ) measured at the same aperture size # ap . The same integrand 
using the RCS background galaxy redshift distribution gives the 
thin dashed line, which is a poor match to galaxy auto-correlation 
function. The integrand for the galaxy-mass cross-correlation func- 
tion, /i3(ui), is indicated by the thick dotted line, which overlaps 
well with the galaxy auto-correlation function. 



h 3 (w;9 ap )= , Pmtcr(w]9 ap ). (16) 

a(w)f K (w) 

The parameter r in Eqn. 15 is the galaxy-mass cross- 
correlation coefficient (e.g., Pen 1998; Dekel & Lahav 1999; 
Somerville et al. 2001), which is a measure of non-linear 
stochastic biasing. In the case of deterministic linear bi- 
asing, r = 1. 

The redshift distributions of the lenses and the sources 
used to measure the galaxy-mass cross-correlation function 
partly overlap. A source in front of a lens contributes no 
signal, and therefore lowers the lensing signal. The overlap 
of the redshift distributions is naturally taken into account 
by Eqn. 15. 

To prove this statement, we first consider a sheet of 
lenses at a distance uif. The "lensing strength" for these 
lenses is given by g(w/), to which only background galaxies 
with lo > ujf contribute. The larger the u)f, the lower both 
the observed and predicted lensing signals will be (and 
they are lower by the same factor). Hence, Eqn. 15 gives 
the correct result for a sheet of lenses at a given redshift. 
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The actual redshift distribution of lenses can be consid- 
ered a superposition of many sheets at different rcdshifts, 
and consequently Eqn. 15 is correct for any combination 
of lens and source redshift distributions. 
The combination of Eqns. 7, 11, and 15 relates the value of 
r to the observed correlation functions and the cosmology, 
and is given by 

\J I dwh! (w;9 ap )^J / dwh 2 (w; 9 ap ) (M ap 7V) 



Jdwh 3 (w,9 ap ) (M a 2 p ) 1/2(^2) 1/2 



= h(0a.p, Mm, ^A) X 



{M ap (9 ap )M(9 ap )) 



(17) 



WM><A4 2 p(M> 

Similar to /i, the value of / 2 depends minimally on 9 ap 
and the assumed power spectrum (Van Waerbeke 1998). 

The integrand h^(w;9 ap ) of the integral in Eqn. 15 is 
also shown in Figure 1, and it matches the results for the 
galaxy auto-correlation function well. As for the bias pa- 
rameter 6, the value of r can depend on scale and redshift. 
However, because of the properties of the aperture statis- 
tics, described in Appendix A, we can measure the scale 
dependence of r directly from Eqn. 17. As for 6, the in- 
ferred value for r is a redshift averaged value for the lenses 
in our sample. 

Our definitions of b and r are chosen such that they 
can be related directly to the observed correlation func- 
tions. However, both b and the cross-correlation coefficient 
r mix non-linear and stochastic effects (Dekel & Lahav 
1999; Somervillc et al. 2001), which we currently cannot 
disentangle with weak lensing. 

Dekel & Lahav (1999) defined a combination of 3 pa- 
rameters (6, 6, and erf,), which separate non- linear and 
stochastic biasing. They parametrize the stochasticity by 
the parameter <T(,, the local biasing scatter. The non-linear 
biasing is characterized by the ratio 6/6, where 6 is the 
slope of the linear regression of 5 g on 5 (i.e., the natural 
generalization of the linear biasing parameter) . 

From lensing we can only measure combinations of these 
parameters. We can, however, relate the observables b 
and r to the the parameters from Dekel & Lahav (1999) 
in the two limiting cases. If we assume that the bias is 
purely non- linear and deterministic (i.e., at = 0), we ob- 
tain 6 = 6, and 6/6 = 1/r. Hence, the inverse of the 
cross-correlation coefficient measures the amount of non- 
linear biasing. Likewise, in the case of linear and stochastic 
biasing (i.e., b — 6) we obtain <t& = b\J\ — r 2 . 

4. DETERMINATION OF THE CORRELATION FUNCTIONS 

A straightforward implementation of the method de- 
scribed above is to tile the survey area with apertures, and 
compute (M ap N){9) and (Af 2 )(9) directly from the data. 
To do so, one can use the estimators (Schneider 1998): 



M a P = 7rfl ap + , and N = ^ U(6i), 

l^i=l Wi i=l 

where Nf and Nb are respectively the number of lens and 
source galaxies found in the aperture of radius 9 ap , N is 
the average number density of lenses, and j t ,i is the ob- 
served tangential shear of the zth background galaxy. The 



weights Wi correspond to the inverse square of the uncer- 
tainty in the shear measurement (Hoekstra et al. 2000). 

However, this procedure has the disadvantage that it 
assumes a contiguous data set; i.e., there are no holes in 
the data. In real data, regions that are contaminated by 
bright stars, bad columns, etc., have to be masked. Al- 
though the masking of the RCS data is not severe, we will 
use a different approach, which is much less sensitive to 
the geometry of the survey. 

Instead we measure the ensemble averaged tangential 
shear as a function of radius around the sample of lenses 
("galaxy-galaxy lensing" signal) and use these results to 
derive (M ap Af). We use the angular two-point correla- 
tion function to estimate (Af 2 ). The mass auto-correlation 
function (M 2 p ) is derived from the observed ellipticity cor- 
relation functions (e.g., Pen et al. 2002; van Waerbeke et 
al. 2002; Hoekstra et al. 2002b). We show below how 
these observed correlation functions, which do not require 
contiguous survey areas, can be easily related to the aper- 
ture mass correlation functions. 

We first consider the angular two-point correlation func- 
tion uj(9), which is related to the power spectrum through 



2tt 



cj(9) = (An g (0)An g (9)) 
2 

' " ' U 3d 



M^m) ! illp Ajhr w ) 3 ° m - (18) 

where Jo (x) is the zeroth order Bessel function of the first 
kind. We can define the effective (projected) power spec- 
trum of the angular correlation function P u (l) as 



dw 



p f (w) 



P: 



I 



(19) 



Ik{w)J ' \fK{w) 

Thus we obtain for the angular correlation function 

u{9) = ^J dl IPuWMW)- (2°) 

For the auto-correlation function of Af we have (see Eqns. 7 
and 9) 



(JV 2 }{9) = 2tt / dl IP U (1) 



12J 4 (19) 



-1 2 



(21) 



Schneider, van Waerbeke & Mcllicr (2002) have shown 
that it is rather straightforward to transform one corre- 
lation function into another. Using the orthogonality of 
Bessel functions, we obtain 



PUi) 



= 2. / 



d$ tfw(tf)J (W), 



(22) 



which we can use to relate (Af 2 )(9) to lo{9). Doing so, we 
obtain 



■d 



(Af 2 )(9)= I d-d - u (0)T+[- 



(23) 



where the function T+(x) is defined as (Schneider ct al. 
2002) 

T+(x)- 576y o _J (xt)[J 4 (t)] 2 . (24) 
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Schneider et al. (2002) found an expression for T + (x) in 
terms of elementary functions 



T + {x) = 



6(2 - 15a; 2 ) 



1 arcsin(x/2) 

7T 



-(120 + 2320a; 2 - 754a; 4 + 132a; 6 - 9a; 8 ), (25) 



IOOtt 

for x < 2. T + (x) vanishes for x > 2. Therefore the integral 
in Eqn. 19 extends only over < < 20. Hence we need 
to measure the angular correlation function out to twice 
the aperture size in order to compute (A/" 2 )(0). 

We now turn to the measurement of the galaxy-mass 
cross-correlation function. In this case the azimuthally av- 
eraged tangential shear around the lens galaxies is a very 
useful statistic (e.g., Fischer et al. 2000). To derive the 
galaxy-mass correlation function we use the relation be- 
tween the average convergence /t(0) inside a circular aper- 
ture of radius and the mean tangential shear along the 
boundary of the aperture 



WW 



1 dR(9) 

2 d\n6' 



(26) 



The galaxy-mass correlation function is defined as 
(An s (O)7t(0)) = (7t)(0), and is related to the power spec- 
trum as (Kaiser 1992; Guzik & Seljak 2001) 



br f dw i^Pr 



4tt 
dl lP 3d 



)Pf(w) 
(w) 



I 



... , ■ : ■ (27) 

Jk{w) ) 

where J 2 (x) is the second order Bessel function of the first 
kind. We assume that the cross power spectrum is related 
to the matter power spectrum by brP 3d (k). As discussed 
above, the assumption of a constant value of b and r does 
not change the interpretation of our measurements. 

We define the effective (projected) power spectrum for 
the galaxy-mass correlation function P 7 (Z) as 



3 fH { 



/ 



dw 9M P/ M p 3d 



1 



which gives 



<7t>(0) 



a{w)f K {w) \fxiw 



±-Jdl lP 7 (l)h(W), 



and 



(AfM ap )(0) = 2tt J dl LP 7 (Z) 



12J 4 (Z0) 



(28) 



(29) 



(30) 



As before we can express P 7 (Z) by an integral over the 
galaxy-mass cross-correlation function, and we obtain 



(AfM ap )(9)= J dd ^(>ytW)F 
where we define the function F(x) as 



(31) 



F{x) = 576 [ 
Jo 



p 2 (xt)[J 4 (t)} 2 . (32) 



We did not find a simple expression for this function, 
but it is well behaved for < x < 2, and vanishes for 
x > 2. Hence, as before, the integral in Eqn. 31 extends 
only over < < 20. 

The mass auto-correlation function (Mf p ) can be ob- 
tained by measuring the ellipticity correlation functions, 
which are given by 



Eij w i^7t,i(xi)-7t,j(xj) 



£tt(0) = 



6r(0) = 



V^b 

Ei,, W i W j 



v^b ' 
Eij WiWj 



(33) 



(34) 



where = |xj — Xj|. j t and 7 r are the tangential and 45 
degree rotated shear in the frame defined by the line con- 
necting the pair of galaxies. For the following, it is more 
useful to consider 

C+(0) = £ tt (0) + M*)> and = fttO?) -M0), (35) 

i.e., the sum and the difference of the two observed cor- 
relation functions. As shown by Crittenden et al. (2002), 
one can derive "E" and "B"-mode correlation functions by 
integrating £+(0) and £-(0) with an appropriate window 
function (see Pen et al. 2002 for an application to the 
VIRMOS-DESCART data). 

The "E" and "B"-mode aperture masses are computed 
from the ellipticity correlation functions using (e.g., Crit- 
tenden et al. 2002; Schneider et al. 2002) 



<M 2 p >(0) 



and 



(M|)(0) = 



dd g 



dd_d_ 



t) 



d 



(36) 



7 ap 



-€_(<?)T_ 



'ap 



(37) 

The expression for T + (x) is given by Eqn. 22, whereas 
T_(x) is given by 



2\ 7 /2 



192 



(38) 



for x < 2, and T_(x) vanishes for x > 2. 

The "B"-mode aperture mass (Mf) provides a quantita- 
tive estimate of the systematics, since gravitational lensing 
only produces an "E"-mode. Residual systematics (e.g., 
imperfect correction for the PSF anisotropy) or intrinsic 
alignments will give rise to a "B"-mode. A similar test 
exists for the galaxy-mass cross-correlation function: in 
the absence of systematics, the average signal around the 
lenses, when the sources are rotated by 45 degrees, should 
vanish. 
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Fig. 2. — The measurements of (W 2 ) (panel a), and (ATM ap ) 
(panel b) as a function of angular scale from the RCS data. The 
largest scale corresponds to an aperture radius of 1 degree. Panel c 
shows (M 2 p ) as a function of angular scale from the VIRMOS- 
DESCART data (largest scale is 45 arcminutes). The error bars 
for (Mf p ) have been increased to account for the unknown correc- 
tion for the observed "B"-mode. For reference (they are not fitted 
to the measurements) we have also plotted model predictions (for 
6=1 and r = 1) for an OCDM cosmology (dotted line; Q m = 0.3, 
Cl A = 0, cr 8 = 0.9, and T = 0.21) and a ACDM cosmology (dashed 
line; H m = 0.3, Q A = 0.7, a 8 = 0.9, and T = 0.21). Note that the 
points at different scales are slightly correlated. 



5. MEASUREMENTS 

In this section we present the measurements of 
the galaxy auto-correlation and the galaxy-mass cross- 
correlation function (M ap (# ap ).A/'(# a ,p)), which were ob- 
tained from the RCS data. The mass auto-correlation 
function (M 2 p (6> ap )) was measured by van Waerbeke et 
al. (2002) from the VIRMOS-DESCART survey. 

As described in §4, we do not measure (A -2 ) directly 
from the data. Instead we measure the angular correla- 
tion function w{9) and use Eqn. 23 to determine (A/" 2 ). To 
measure ui(6), we use the well known estimator (Landy & 
Szalay 1993) 



u{0) = 



DD - 2DR + RR 
RR ' 



(39) 



where DD, DR, and RR are pair counts in bins of 9 + 68 of 
the data-data, data-random, and random-random points 
respectively. For each patch, we create 24 mock catalogs 
by placing the lens galaxies at random positions in the un- 
masked regions of the data. We note our sample of lenses 
is complete within the unmasked regions: the galaxies are 
sufficiently faint that they are not saturated, and they are 
several magnitudes brighter than the detection limit. Fur- 
thermore, the uncertainties in the photometry are small 
(< 3%, Gladders et al. in preparation), and do not affect 
the measurement of the angular correlation function of the 
lens galaxies. 



In order to determine the true angular correlation func- 
tion one needs to apply an integral constraint correction 
to the observed correlation function. The determination 
of this correction is not trivial, and can introduce a signif- 
icant systematic uncertainty. However, since the integral 
J dx xT + (x) vanishes, (A/ -2 ) has the nice property that it is 
independent of the integral constraint correction. In addi- 
tion, the measurements on different scales are only mildly 
correlated. 

Although we use the observed angular correlation func- 
tion as an intermediate step in our analysis, it is use- 
ful to compare the results to previous work. The ob- 
served angular correlation function is well approximated 
by a power law with slope -0.7, and an amplitude w(l') = 
0.115±0.005. Unfortunately, a direct comparison with the 
literature is difficult, because of the different sample selec- 
tion (magnitude limits, filters). We can, however, make 
a crude comparison with the results from Postman et al. 
(1998), who measured the angular correlation function in 
the /-band. The mean i?c-band magnitude of our sample 
is ~ 20.5, which corresponds to a mean /-band magnitude 
of ~ 19.8 (based on the colors of galaxies observed in the 
CNOC2 survey). Postman et al. (1998) list a value of 
uj(V) = 0.136 ± 0.021 for galaxies with a median magni- 
tude of I = 19.6, which is in reasonable agreement with 
our result. 

We measured (A -2 ) for each of the 10 RCS patches, and 
used the scatter in these measurements to estimate the 
error bars on the galaxy auto-correlation function. The 
errors therefore also include cosmic variance. The results, 
measured out to 1 degree, are presented in Figure 2a. 
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Fig. 3.— (a) The observed ratio (M ap A/"> / '((•Af 2 >( M fp» 1/2 as a 
function of aperture size, (b) The ratio of the mass auto-correlation 
(M 2 p ) and the galaxy auto-correlation (Af 2 ) as a function of aper- 
ture size. The measurements at different scale are slightly corre- 
lated, and the error bars correspond to the 68% confidence inter- 
vals. The upper axis indicates the effective physical scale probed 
by the compensated filter U{<f>) at the median redshift of the lenses 
(z = 0.35). The dotted lines are the predictions of an OCDM model, 
whereas the dashed lines correspond to the ACDM model for b = 1, 
and r = 1. In both cases, if b, and r are constant with scale, we 
expect to observe ratios that are virtually constant. The observed 
ratios show significant variation with scale, thus implying that both 
b and r vary. 
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Fig. 4.— (a) The observed ratio {M x M)/({M 2 )(M 2 p )) 1 / 2 as a 
function of aperture size. This corresponds to the results when the 
phase of the shear is increased by n/2. If the signal presented in Fig- 
ure 3a is caused by lensing it should vanish here, as it does, (b) The 
ratio of the "B"-mode auto-correlation (M 2 .) and the galaxy auto- 
correlation (A/" 2 ) as a function of aperture size. The results suggest 
a residual "B"-mode, which is attributed to systematics (note that 
the points at different scales are somewhat correlated). The ampli- 
tude of this signal is added to the error estimates in Figure 3b as 
a conservative limit. We note that the amplitude of the "B"-mode 
signal is small compared to the lensing ("E"-mode) signal. 



The observed ratio (M ap A/')/((A/' 2 )(Mf p )) 1 / 2 as a func- 
tion of aperture size is presented in Figure 3a. Figure 3b 
gives the ratio of (Mf p ) and (TV 2 ). The errorbars on 
the ratios correspond to the 68% confidence limits, and 
have been determined from a Monte Carlo simulation us- 
ing the uncertainties in the measurements of the observed 
correlation functions (which were assumed to be Gaus- 
sian). For reference, we have also indicated the effective 
physical scale (~ FWHM of the filter function) probed by 
the compensated filter U((f>) corresponding to a redshift of 
z = 0.35. 

For reference, Figure 3 also shows the ratios for an 
OCDM and a ACDM model in the case 6=1 and r = 1 
(i.e., we have plotted x /i, and / 2 ). These model values 
are virtually constant with scale if b and r are constant. 
Also note that the value for the cross-correlation is almost 
the same for both cosmologies, whereas the ratio of the 
auto-correlation functions differs by almost a factor 2. 

To examine possible systematic effects, we also com- 
puted the results when the sources are rotated by 45°. 
This signal should vanish in the case of lensing. The results 
of this powerful test for the cross-correlation are presented 
in Figure 4, and are indeed consistent with no signal. The 
results for the auto-correlation do show some residual sys- 
tematics. The amplitude of this signal is added to the 
error estimates in Figure 3b as a conservative limit (van 
Waerbeke et al. 2002). We note that the amplitude of the 
"B"-mode signal is small compared to the lensing ("E"- 
mode) signal. Based on these results we conclude that the 
accuracy of our measurements (in particular the galaxy- 
mass cross-correlation) is not limited by systematics. 



Figure 2b shows the cross-correlation function (M ap Af) 
as measured from the RCS. To derive the cross-correlation 
function we measured the azimuthally averaged tangential 
shear around the lens galaxies in bins of 1 arcsecond. We 
then used Eqn. 31 to relate the observed tangential shear 
profile to (M ap Af). As before, we measured the cross- 
correlation for the 10 RCS patches, and used the scatter 
in the measurements to estimate the errorbars. 

The mass auto-correlation function (M 2 p ), presented in 
Figure 2c, was taken from van Waerbeke et al. (2002). ft 
was derived from the observed ellipticity correlation func- 
tions (see Eqn. 36). The largest scale measurement of 
(Mf p ) from the V1RMOS-DESCART survey is 45 arcmin- 
utes. The error bars in Figure 2c have been increased to 
account for the unknown correction for the observed "B"- 
mode (van Waerbeke et al. 2002). 

The shapes and amplitudes of the correlation functions 
presented in Figure 2 depend on the power spectrum. As 
discussed in §3, we can remove the dependence on the 
power spectrum by taking ratios of the correlation func- 
tions. Since the various correlation functions probe the 
power spectrum at the same redshift (see Fig. f), we can 
take the ratios of measurements at the same angular scale. 
If the surveys were not well matched, we would have had 
to compare the measurements at different angular scales, 
in order to ensure that we probe the power spectrum at 
the same physical scale. In addition, one has to take into 
account the evolution of the power spectrum in such a 
situation. 



6. DISCUSSION 

The results presented in Figure 3 suggest significant 
variation of both b and r with scale. We convert the ob- 
served ratios into estimates for the bias parameters for the 
currently favored cosmological model (£l m = 0.3, Ha = 
0.7), using Eqns. 14 and 17. 

Figure 5a shows the inferred value of the galaxy-mass 
cross-correlation coefficient r as a function of scale. Be- 
cause of the small value of (M 2 p ) at I arcminute, the un- 
certainties in both b and r are very large, and we have 
omitted this point. 

We find that on small scales the measurements are con- 
sistent with a value r <~ 1. This good correlation between 
mass and light on small scales (i.e., around galaxies) indi- 
cates that (luminous) galaxies are surrounded by massive 
halos. We note that with our definition r can be larger 
than 1. On the largest scales the results are consistent 
with r = 1. On scales ~ l^sn 1 Mpc r is significantly lower 
than unity, with a minimum value of r = 0.57^q'q^ (68% 
confidence). 

Figure 5b shows that the scale dependence of b is very 
similar to that of r (also see Figure 6). We find that b 
is smaller than unity on scales <~ 1 — 2/i^q 1 Mpc, with a 

minimum value of b = 0.7lt£u5 (at 1.5h^ Mpc; 68% 
confidence) . Hence, the dark matter is more strongly clus- 
tered than the galaxies. The variation of b with scale is 
significant, but a better determination of (M 2 p ) is needed 
to study the scale dependence in more detail. 
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Fig. 5. — (a) The measured value of the galaxy-mass cross correla- 
tion coefficient r as a function of scale for the ACDM cosmology, (b) 
The bias parameter b as a function of scale. The upper axis indicates 
the effective physical scale probed by the compensated filter U(<j>) 
at the median redshift of the lenses (z = 0.35). The errorbars corre- 
spond to the 68% confidence intervals. Note that the measurements 
at different scales arc slightly correlated. 

Our result is in good qualitative agreement with the 
findings of Jenkins et al. (1998) who combined the ob- 
served APM two-point correlation function and the matter 
correlation function derived from numerical simulations. 
We note that a direct comparison cannot be made with 
other measurements (such as those from Jenkins et al. 
1998), because the bias properties depend on galaxy type 
and redshift. For instance, observations show that the 
amplitude of the galaxy correlation function depends on 
luminosity (e.g., Benoist et al. 1996; Norberg et al. 2001). 

Our results are derived for a specific subset of galaxies, 
and are redshift averaged values (the lenses have a large 
range in redshift, with intrinsically brighter galaxies at 
higher redshifts). This is clearly demonstrated by Eqn A7. 
In the simple case that the average bias parameters change 
approximately linearly with redshift, our results can be in- 
terpreted as a measurement of the biasing properties of L* 
galaxies at redshift ~ 0.35. However, the redshift depen- 
dence of the bias in magnitude limited samples is usually 
more complex. 

Hoekstra et al. (2001a) measured the ratio b/r on angu- 
lar scales out to 12.5 arcminutes, and found that the results 
were consistent with a constant value. With the additional 
data (which probe larger scales and give smaller errorbars) 
we find evidence for a small trend with scale. The measure- 
ments of the current RCS data extend out to an angular 
scale of 1 degree (which corresponds to a physical scale of 
12.5/i^q 1 Mpc), and are presented in Figure 6a. Figure 6 
shows that even on a 1 degree scale, the systematics are 
much smaller than the signal. Hence, the value of b/r can 

3 http://www.cfht.hawaii.edu/Science/CFHLS/ 



be determined accurately. The new data yield an average 
value of b/r = 1.090 ± 0.035. Limiting the measurements 
to the angular scales studied in Hoekstra et al. (2001a), 
we find an average value of b/r = 1.05 ± 0.04, in excellent 
agreement with the previous estimate of b/r — 1.05lg '\q. 

Our measurements should be compared to models of 
galaxy formation. The two commonly used approaches 
are hydrodynamic simulations (e.g., Blanton et al. 2000, 
Yoshikawa et al. 2001) or semi-analytical models (e.g., 
Kauffmann et al. 1999a, b; Somerville et al. 2001, Guzik 
& Seljak 2001). These studies suggest values of r close 
to unity. Unfortunately, the mass resolution in these sim- 
ulations, such as the GIF results (e.g., Kauffmann et al. 
1999a, b), is too poor: they only resolve massive, luminous 
galaxies. Our sample of lenses contains many lower lumi- 
nosity systems, which seriously hampers the comparison 
of our results with predictions. 

As the value of r is intimately linked to the details 
of galaxy formation, a careful comparison of the models 
with the weak lensing measurements provides unique con- 
straints. In addition, planned large weak lensing surveys, 
such as the CFHT Legacy Survey 3 , will significantly im- 
prove the accuracy of the measurements, allowing larger 
scales to be probed, and enabling us to select galaxies 
based on the colours or luminosities (using photometric 
redshifts) . Therefore the prospects of constraining biasing 
parameters from weak lensing are excellent. 
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Fig. 6. — (a) The value of b/r as a function of angular scale, 
under the assumption Q m = 0.3 and Q/± = 0.7. These results are 
based on RCS data only, which allows this ratio to measured out to 
1 degree (which corresponds to a physical scale of 12.5/i^ Mpc). 
Note that the points are slightly correlated. The value of b/r is 
almost constant over the range probed here. For this cosmology we 
find an average value of b/r = 1.090 ± 0.035, in excellent agreement 
with the result from Hoekstra et al. (2001a). (b) Measurement of 
(M^]\[) / (N 2 ) (the signal when the phase of the shear is increased 
by tt/2), which should vanish if the results in (a) are caused by 
lensing. 
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APPENDIX 

SCALE AND REDSHIFT DEPENDENT BIAS PARAMETERS 

As has been demonstrated in this paper, the bias relation is more complicated than the simple case of linear deterministic 
biasing, and the inferred values of b and r vary with scale. In addition, these parameters are expected to vary with redshift. 

In §3 we treated both b and r as constants, which is not warranted by the data. In this appendix, however, we show 
that as long as b and r vary slowly with scale, we can still infer the bias parameters as a function of scale directly This 
procedure works, because the aperture mass statistic is effectively a pass-band filter. 
If we assume that S g (k) = b(k)6(k), we have to replace Eqn. 5 with 



An g (0)= J dw pf(w)b(f K (w)6;w)5(f K (w)6;w), (Al) 
which changes the galaxy auto-correlation function to 



where 



f PA W ) - 

(AA 2 (0ap)) = 27r / dw-ty-P mei (w-,e ap ), (A2) 
J Jk( w ) 

PfiltcrKM =jdllb 2 (jl-^-w^ P 3d J 2 (W ap ). (A3) 

The filter J 2 (t]) is strongly peaked, and this motivates the approximation of J 2 by a Dirac delta function (Bartelmann & 
Schneider 1999) 

J 2 (V) ~ lrr J (r] - 693tt/512) w 1.43 x lO" 2 ^?? - 4.25). (A4) 
ll557r d 

Consequently Pfiitor can be approximated as 
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^filter (W, ap ) W fr 



4.25 



fl ap /x( w ) ' 



; mi 



-Pfilter(w;flap), 



where Pmter is given by Eqn. 9. For the ratio of the galaxy and matter auto-correlation functions we can write 



(M 



2 ) 



He 







/ dwh 2 (w; 6» ap ) 



/ dwh^w; flap )6 2 (4.25/ (fl ap /^ (tu) ); w) 



(A5) 



(A6) 



or 



/ dwM«;;0a P )& 2 (4.25/(flap/KH); W ) _ 9 /JT, 



/ dwhi(w;6 ap ) 



- I — I n; 



Jrf W fe 2 ( W ;fla P ) (A/" 2 ) 
7^ 1 ( U ;;flap)(M a y 



2 (' /V ' 2 ) 
/l x X ^2 J (^ a p)- 



(A7) 



The functions hi, and hi are given by Eqns. 8 and 12. The actual calculations show that /i does not depend on the power 
spectrum (e.g., van Waerbeke 1998). The left hand side of Eqn. A7 shows that we measure the bias parameter weighted 
by the function h\. 

A similar result can be obtained for the galaxy-mass cross-correlation coefficient r 



jdwh 3 {w;6^)r{4.2h/{e^f K {w));w) _ yj 1 dwh i K M yf f dwh 2 (w; fl ap ) 



(M ap AT) 



J dwh 3 {w;6 a -p) 



h x 



/ dwh 3 (w;O ap ) 
(M ap (fl ap )AA(fl ap )) 

^(^(fl a p))(M a 2 p (flap))' 



(A8) 



The redshift distribution of the lenses Pf(w) is narrow (see Fig. 1). Hence, we measure the bias parameters over a small 
range in redshift. Because Jk{w) varies somewhat over this interval, we average values of b and r on different physical 
scales. However, if 6 and r vary slowly with scale and redshift, we can ignore the variation of these parameters over this 
range in scale. Hence, Eqns. 12 and 14 can be used to obtain a direct measurements of the bias parameters as a function 
of scale. 



